Probing a topological quantum critical point in semiconductor-superconductor 

heterostructures 



O 
(N 



u 



I 

o 
o 



O 

in 
o 




Sumanta Tewari-'^, J. D. Sau^, V. W. Scarola'^, Chuanwei Zhang**, and S. Das Sarma^ 
^Department of Physics and Astronomy, Clemson University, Clemson, SC 29634 USA 
^ Condensed Matter Theory Center and Joint Quantum Institute, Department of Physics, 
University of Maryland, College Park, Maryland, 20742-4111, USA 
'^Department of Physics, Virginia Tech, Blacksburg, VA 24061, USA 
'^Department of Physics and Astronomy, Washington State University, Pullman, WA 99164 USA 

Quantum ground states on the non-trivial side of a topological quantum critical point (TQCP) 
have unique properties that make them attractive candidates for quantum information applications. 
A recent example is provided by s-wave superconductivity on a semiconductor platform, which is 
tuned through a TQCP to a topological superconducting (TS) state by an external Zeeman field. 
Despite many attractive features of TS states, TQCPs themselves do not break any symmetries, 
making it impossible to distinguish the TS state from a regular superconductor in conventional bulk 
measurements. Here we show that for the semiconductor TQCP this problem can be overcome by 
tracking suitable bulk transport properties across the topological quantum critical regime itself. The 
universal low-energy effective theory and the scaling form of the relevant susceptibilities also provide 
a useful theoretical framework in which to understand the topological transitions in semiconductor 
heterostructures. Based on our theory, specific bulk measurements are proposed here in order to 
characterize the novel TQCP in semiconductor heterostructures. 



Introduction: Quantum critical points (QCP) sepa- 
rate two many-body quantum ground states distinguish- 
able by a macroscopic order parameter M (see Fig. [T^) 
Q. In Fig.UJi the solid curve denotes a true phase tran- 
sition line while the dashed curve, given by fc^T ~ Eq 
where Eq is the zero-temperature energy gap, represents 
only a crossover. The two curves meet at the QCP at a 
specific value of the tuning parameter g = gc, straddling 
a finite regime in the (T — (7) plane usually called the 
quantum critical (QC) regime 0] . Quite interestingly, at 
the QC regime the effects of the zero-T quantum phe- 
nomena and the associated QCP are manifest even at 
finite temperatures. This so-called 'quantum fan' region, 
where quantum criticality manifests far from g = gc &t 
finite temperatures, enables experimental studies of QCP 
which are strictly speaking T = phase transitions tuned 
by the parameter g. 

Topological quantum critical points also separate, 
based on distinct topological properties, two macroscopic 
quantum ground states, although the states in question 
now have exactly the same symmetries and thus cannot 
be distinguished by any local order parameter or bulk 
measurements (Fig. [TJj) Q. Frequently, the quantum 
state on the topological side of a TQCP can be distin- 
guished by certain non-trivial statistical properties of its 
excitations 043 j ^^^^ ^ novel ground state quan- 
tum degeneracy which is not associated with any sym- 
metry in the underlying Hamiltonian Q. An example 
is provided by an electron- or hole-doped semiconduc- 
tor thin film or nanowire with s-wave superconductivity 
tuned through a TQCP by an externally applied Zeeman 
splitting P. This system has recently been studied exten- 
sively after it was pointed out by Sau et al. |i8] that for 
P greater than a critical value Pc this system supports 



FIG. 1. (a) Phase diagram associated with a conventional 
QCP in the (T — (;)-plane, where g is the tuning parameter 
in the Hamiltonian. Solid curves denote true phase transi- 
tions, while dashed curves denote only a crossover, (b) Finite- 
temperature phase diagram associated with the TQCP in a 
spin-orbit coupled semiconductor. The tuning parameter F 
represents a suitably directed Zeeman splitting. The super- 
conducting pair potential A is perfectly continuous and non- 
zero at, and on either side of, the TQCP (T = 0,F = Fc). 
Consequently, all lines on this diagram are crossover lines de- 
fined by fcflT — Eo where Eo is the zero-temperature energy 
gap. 



novel non-Abelian topological states [9|-|20[. For P > Pc, 
defects in the (proximity-induced) s-wave pair potential 
A can support localized topological zero-energy excita- 
tions called Majorana fermions. Majorana fermions, with 
second-quantized operators 7 satisfying 7^ = 7, follow 
non-Abelian exchange statistics under pair- wise exchange 
of the coordinates [4|-l7[. Majorana fermions have been 
predicted to be useful for building a topological quantum 
computer which is intrinsically fault-tolerant to all local 
environmental decoherence 0, 0] ■ 

The existence of Majorana fermions in the defects of 
A notwithstanding, A itself remains perfectly continuous 
and non-zero on both sides of the TQCP in a semicon- 
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ductor. This leads to there being no quahtative differ- 
ence between the two states in conventional bulk mea- 
surements. Thus, it appears that simple transport quan- 
tities such as resistance are unable to demonstrate the 
emergence of the TS state. In this paper we propose a 
very specific (and experimentally simple) scheme for the 
direct observation of the TQCP in bulk measurements 
provided such measurements access the so-called topo- 
logical quantum critical regime, which occurs before the 
system settles into the TS state at large F > Fc. We 
believe that our work brings topological quantum phase 
transitions (TQPT) in semiconductors explicitly into the 
mainstream of quantum critical phenomena, something 
that was only implicit in the extensive existing literature 
on this class of TQPTs in semiconductors. 

Hamiltonian, TQCP, and phase diagram at fi- 
nite temperatures: The semiconductor {e.g., InAs) 
system mentioned above is mathematically described 
by the following Bardcen-Cooper-Schrieffer (BCS)-type 
Hamiltonian; 

H = {rjk^ - /^)t^ + FS • (t + ^(fe X (t) • ir^ + Ar^^ (1) 

where S' is a suitably chosen direction of the applied Zee- 
man spin splitting given by F = 1/2(7 /x^i? with g the ef- 
fective Lande g-factor, B the applied magnetic field and 
fjLB the Bohr magneton. H is written in terms of the 
4-component Nambu spinor (u-|-(r), U4.(r), W4,(r), —v^{r)), 
and the Pauli matrices <Jx,y,z,'''x^y^z act on the spin and 
particle-hole spaces, respectively. H can describe a 2D 
system when fc = {kx,ky) is a 2D vector with S = z, 
while a ID structure is described by choosing k = 
with S = X. Here, rj = 1/m* with to* the effective 
mass of the charge-carriers, fj, is the chemical potential 
measured from the bottom of the top-most confinement 
induced band, the Zeeman splitting F breaks the time re- 
versal symmetry, a is the Rashba spin-orbit coupling con- 
stant, and A is an s-wave superconducting pair-potential 
proximity induced in the semiconductor from an adjacent 
superconductor {e.g., Al). 

The Hamiltonian in Eq. ([T]) has recently been stud- 
ied extensively [sl-fioj. A TQCP exists in this system 
as the tuning parameter F is varied through the crit- 
ical value F = Fc = a/A^ -f /i^ where the quantity 
Co = (A^ + - F^) changes sign. For Co > 0, the (low- 
F) state is an ordinary, non-topological superconductor 
(NTS) with only perturbative effects from the Zeeman 
and spin-orbit couplings. For Co < 0, however, the (high- 
F) state has non-perturbative effects from a, and can sup- 
port zero-energy Majorana fermion excitations localized 
at the defects of the pair-potential A Q . The parameters 
for the TQCP are in the experimentally achievable range 
because, for a typical semiconductor wire say InAs or 
InSb, because of a large effective g ^ 15 — 50 a moderate 
magnetic field _B ~ 0.5 T corresponds to a Zeeman split- 
ting F '-^ 2 — 8 K. Noting that /i in Eq. ([T]) corresponds to 



the chemical potential measured relative to the bottom of 
the top-most confinement band and A '--^ 1 — 10 K for an 
ordinary s-wave superconductor, a moderate B ^ 0.5 T 
should be sufficient to induce the topological phase tran- 
sition in the semiconductor. For numerical calculations 
in this paper we have assumed fj, — (Fermi surface at 
the bottom of the top-most band), A = 0.5 meV, a = 0.3 
meV, so that Fc ~ 0.5 meV. The existence of Majorana 
fermions at defects, and also the fact that the high-F 
state is isomorphic to a spin- less px + ipy superconductor 
lol 21, 22 1, make the high-F state a topological super- 
conductor. Interestingly, A remains non-zero and con- 
tinuous across the TQCP [H [l^, so the NTS and TS 



states break exactly the same symmetries, namely, gauge 
and time-reversal. As a result, no macroscopic local or- 
der parameter can differentiate between the NTS and TS 
states, and they cannot be distinguished by any known 
bulk measurements. 

For our present purposes note that the topological crit- 
ical point Fc is marked by the single-particle minimum 
excitation gap Eq vanishing as a function of the Zeeman 
splitting. This can be seen by diagonalizing the Hamil- 
tonian in Eq. ([T]) to obtain the lower-branch of the quasi- 
particle excitation spectrum. 
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where e — rjk^ ~ /i and r| = F^ 



For F near Fc, 



the minimum of _Bfc is at fc = 0, and setting fc = in 
Eq. ^ we find the minimum quasiparticle gap Eq given 
by, 



Eq = \t- Va^+^I 



(3) 



Eq vanishes exactly at Co = 0, which marks the TQCP 
separating the NTS and TS states. Note that Eq is finite 
and positive for both Co > (NTS state) and Co < 
(TS state). 

It is important to note that the system exactly at the 
zero-temperature TQCP (T = 0, F = Fc) can be thought 
of as an s-wave superconductor (A is finite and con- 
tinuous at F = Fc). This bulk s-wave superconductor, 
however, coexists with nodal fermions at fc = 0. Since 
on both sides of F = Fc the ground states are fully 
gapped, we can construct the {T — F) phase diagram 
for this TQCP (Fig. [T]d) by drawing two crossover curves 
marking ksT ^ ^o(F) on both sides of Fc. Note that 
since A is non-zero everywhere on the phase diagram, 
the finite-temperature crossover curves can only be justi- 
fied, as we will show below (see the discussion following 
Eq. [T2)) , on the basis of some measurable quantities show- 
ing a pronounced change across these curves [2^. Such 
identification of crossover curves opens a finite regime in 
the phase diagram which we can associate with the quan- 
tum critical regime of this Zeeman-tuned TQCP. The QC 
regime (and the crossover curves) can also be understood 



2J| as the regime in which the system, with increas- 
ing length scales, encounters the thermal length scale 
(/3 — {kBT)~^) before it encounters the zero-temperature 
correlation length ^ (/? < ^ ^ (F — Fc)"^), which diverges 
as F approaches Fc. 

Bulk measurement (AC response): We now ask 
what sort of bulk measurements can access the nodal 
fcrmion spectrum and reveal the underlying TQCP at 
T = 0. At first glance it may appear that a simple DC re- 
sistance measurement may suffice, because, as a function 
of increasing F at low T, one should first see negligible 
resistance (NTS state with a gapped spectrum) followed 
by nonzero resistance (QC regime with nodal fermions) 
and finally again negligible resistance (TS state with re- 
entrant gapped spectrum). This idea in practice would 
not work, however, because in the presence of a DC volt- 
age, the superconducting condensate will short the cur- 
rent out even in the QC regime, thus producing negligible 
resistance everywhere in the phase diagram with increas- 
ing F. 

We consider an alternative route to identifying the 
TQCP: an AC measurement. For the sake of definite- 
ness, we consider below (Fig. [5]) the AC conductivity 
across a ID nanowire contacted by s-wave superconduct- 
ing leads which produces the proximity effect. For AC 
conductivity measurements, depending on the frequency 
of the applied voltage, there will be significant excitation 
of the nodal fermions only when the energy scales cor- 
responding to the frequency or the temperature become 
of the order of the single particle gap Eq. The excited 
nodal fermions contribute to the dissipative (real) part of 
the complex AC conductivity, while the superconducting 
condensate contributes to the inductive (imaginary) part. 
Thus, as a function of increasing F starting from deep in 
the gapped NTS state (F = 0), one should first observe 
negligible, followed by non-zero, and then re-entrant neg- 
ligible dissipative response in AC conductivities. Such a 
behavior as a function of F marks a cross-over across the 
dashed curves in Fig. [Hd and in turn reveals the under- 
lying zero-temperature TQCP. 

For quasi-2D and 3D systems such a dissipative AC re- 
sponse of the conductivity due to nodal fermions can be 
tracked by electromagnetic absorption experiments [25| . 
In the present case of a ID semiconducting nanowire 
proximity-coupled to bulk s-wave superconductors, such 
absorption experiments can be difficult, and another 
quantity is needed which nevertheless is still given by the 
real part of the AC conductivity suitably defined. Below 
we show that the AC Josephson impedance serves this 
purpose and it can be easily measured across a ID semi- 
conducting nanowire contacted by s-wave superconduct- 
ing leads which also produce the proximity effect (Fig.[2|). 

AC Josephson impedance: The simplest AC re- 
sponse function that can be measured in a ID nanowire in 
proximity contact with superconductivity is the Joseph- 
son impedance. The linear response of the measured 
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FIG. 2. Nanowire geometry for identifying the topological 
phase transition with increasing Zeeman splitting F. The 
nanowire (shown in green) is contacted by two superconduct- 
ing leads (blue) which produce the proximity effect. The 
leads are placed with a finite potential difi'erence V and the 
current / is measured. All quantities V(ci'), iiy) and '^(lS) 
are frequency-dependent. The josephson phase is $(tj) = 
-iV(cj)/a;. 



current lit) to a sinusoidal voltage Vit) — y(w)e'"* is 
written as l{t) = {x^iyi) + iX\iyi))V{^) where (x2(w) -I- 
iYi(a;))~^ is the Josephson impedance of the junction 
[26| in Fig. [2] In the geometry of Fig. [2l when X2 = 0, 
the junction behaves like a conventional non-dissipative 
Josephson junction. However, when xi becomes non- 
zero, the voltage and current cease to be completely 
orthogonal to each other and a finite amount of power 
(J dtl {t)V (t)) is dissipated in the junction. For the sake 
of brevity, below we will refer to X2 simply as the dissi- 
pative susceptibility keeping in mind that in reality this 
is the real part of the inverse impedance function. 

We assume that the time-dependent voltages at the su- 
perconducting leads are given by VL{x,t) — ~V{t)/2 on 
the left lead and Vii{x, t) = V{t)/2 on the right lead. Cor- 
respondingly, the superconducting phases on the right 
and the left leads are (/)_R(f) = —(j>L{t) = $(i)/2, where 
^(t) is the total time-dependent phase difference between 
the leads. As is well known, V{t) and the phase variation 
<I>(t) are related by the Josephson relation V{t) = ^{t) /2. 

The time-dependent BCS Hamiltonian describing the 
nanowire is given by 



Hi — I dxt/j^—d'^ + V{x, i) — /i + Ta-x — iaaydx]ip 

d2;[A(a;)e*'^(^'*VtV'I + h.c] (4) 



where the pairing field A.{x) is proximity-induced and 
therefore is assumed to be non-zero only in the parts 
of the nanowire in direct contact with the supercon- 
ducting leads. ^/;^ — ipl{x,t) creates a fermion in spin 
state s =ti i s-iid spin-indices are implicitly summed over. 
V{x, t) is the voltage difference across the wire. By apply- 
ing the gauge transformation ■ips{x,t) e^^^^'*^^ps{x,t), 
V{x, t) -> V{x, t)+dtA{x, t), A{x, t) -> A{x, t)+dxA{x, t) 
and (t>{x^t) (l>{x,t) + 2A{x,t) to the Hamiltonian in 
Eq. (HI), and by choosing A{x,t) = —4>{x,t)/2, Hi be- 
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+ aay{-idx - A{x,t))]^/j + A(a;)['(/'j?/'I + h.c] (5) 

where A{x,t) — —dx(t>{x,t)/2. Since the voltage dif- 
ference V{x,t) drops smoothly across the junction with 
width W, we choose a phase dependence (l){x,t) — 
^{t){2erf{x/W) — l)/2 so that increases from 

-$(t)/2 for X < to $(t)/2 at X > VF. 

Because of the spin-orbit coupling term in Eq. (j4]), the 
current operator J takes the following modified form: 

J{x, t) = -i {il^'^dxip - dxip'^'ip) - aip^'a^ip. (6) 

Hi can then be rewritten as, 

Hi{t)^Hio+ dxJ{x,t)A{x,t), (7) 



(to linear order in A) where Hiq is Hi at A — 0. 

The conductance of the wire is calculated from the 
current J{xi,t) at a position xi > W (outside the junc- 
tion) in response to a perturbation / dxJ{x,t)A{x,t). In 
the limit of a small junction (W — > 0), we can approx- 
imate the perturbation as $(i)J(0,t)/2 and the mea- 
sured current as J{0+,t). Therefore, choosing a time- 
dependent voltage V{t) = F(cij)e*'^* corresponding to a 
phase $(i) = V {uj)e'^'^* / (iuj) and using the fluctuation- 
dissipation theorem, the real (dissipative) part of the cur- 
rent response function, X2{^) = Re[5J{0+,uj)/SV{uj)] is 
given by 

X2(c^) = -/m[/ die-'-*([J(0,0),J(0,t)])]. (8) 
^ Jo 

We use Eq. ^ to calculate the dissipative susceptibility 

X2- 

Substituting the operators J(0, t) from Eq. ^ and cal- 
culating the relevant matrix elements we obtain, 



dkdk' 



|(n,fc|TO,fc')r(fc + fc') 



[f{En,k) — f{Em,k')] S{E,n^k' — W — i?„,fe) 



(9) 



Here, f{E) = [e^^^ + 1)^^ is the Fermi occupation func- 
tion, and En^k and |n, k) are the Bogoliubov-de Gennes 
(BdG) eigenvalues and eigenstates of i?i,o, respectively. 
This expression is manifestly a bulk property which is 
real and positive only for positive frequencies. 

The dissipative response X2 vanishes sufficiently far 
from the transition, both in the TS and NTS states, for 
uj smaller than the gap Eq . To see this note that the inte- 
grand in Eq. ([9]) is non-zero only when Em,k' — En^k — 
and Em,k' is empty together with En^k being filled. At 
T = 0, this can only happen when the single-particle gap 
Eq is less than uj. The behavior of X2 as a function of 
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FIG. 3. Dissipative susceptibility X2 (Go = 2e^/h) as a 
function of frequency u. If experimentally u « T is sat- 
isfied, X2 should be plotted with the Zeeman splitting, see 
Fig. 13] and the discussion following it. xa increases at lower- 
frequencies as the Zeeman splitting F is tuned through the 
TQCP (F = Fc = 0.5 meV). For F away from F^, in both the 
NTS and TS states the low- frequency {uu ^ Eo{T)) dissipa- 
tive response is negligible. As the single-particle gap -Eo(F) 
closes at fc = 0, the threshold frequency for the dissipative 
susceptibility approaches zero at the TQCP. Parameters are 
A = 0.5 meV, a = 0.3 meV, /i = and T = 0.1 meV cor- 
responding to 1 K. For reference we have used 1 meV = 250 
GHz. 



UJ for fixed values of T is shown in Fig. [31 As is clear 
from this figure, both the NTS and the TS states are 
non-dissipative for frequencies smaller than a threshold 
set by Eq{T). As F is tuned towards Fc from either side, 
Eo{T) decreases and the threshold value of w for the onset 
of dissipation decreases to zero at the TQCP. Therefore, 
in the experiment suggested in Fig. [^J at sufficiently low 
T,oj <^ A and away from the TQCP, one would expect 
the measured current and the applied voltage to be out of 
phase by 7r/2 so that the power dissipated is zero. As one 
approaches the TQCP, at some value of F, w will surpass 
Eo{T) and a component of the current will become in- 
phase with the applied voltage. This will lead to a finite 
power dissipation, signaling the vicinity of the underly- 
ing TQCP. The succession of behavior with increasing 
F - non-dissipative, followed by dissipative, and then re- 
entrant non-dissipative response - of the Josephson cur- 
rent versus voltage is a clear signature of the underlying 
zero-temperature TQCP. The re-entrant non-dissipative 
response for F > Fc is also a clear signal of re-entrant 
high- Zeeman- field superconductivity, which can only be 
topological in nature [20| . 

Scaling of dissipative susceptibility: For analyti- 
cal calculations of the scaling functions we first need to 
derive the low-energy effective theory valid in the vicin- 
ity of the TQCP. To do this we recall that the TQCP 
is given by the minimum excitation gap Eq vanishing 
as a function of F. Near the transition, only one pair 
of eigenstates of the BCS Hamiltonian in Eq. ([1]) van- 
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ishes near fc = 0. Therefore, near the transition, for 
the low-energy effective theory we can 'integrate out' 
the other pair of Bogoliubov eigenstates and focus only 
on the lowest pair that vanishes at r^. The pair of 
eigenstates n — 1,2, whose energies vanish linearly near 
/c ~ 0, form a pair of chiral Majorana fermion operators 
7„(x) = / dke'''^Y.s[^n,s{k)rs{k) + VnAk)i^s{k)]. Here, 
{un,■\{k),Un,i{k),Vn,i{k),-Vn,^;{k)) are the BdG eigen- 
states n = 1,2 with eigenvalues ±Ek (see Eq. [2]) for 
r = Fc. The low-energy effective action valid near the 
TQCP can then be written in terms of the Dirac fermions 
^^{x) = 7i(x) + ^72 (a;) as, 



dT / dx[^'ta^* + iw(*1'a^^t + /j c) + J*^*], 



(10) 

where 8 — (F — Fc) is the gap which takes the system 
away from the phase transition and v = ab? -I- A^) 
is a velocity determined by the spin-orbit coupling con- 
stant a. It follows that the dynamic critical exponent 2, 
which relates the spatial and temporal correlation lengths 
^ and by ^ = is 1. Since the energy gap b vanishes 
linearly with F, ^ also diverges linearly with the 
Zeeman coupling, and therefore the mean field correla- 
tion length exponent v = \. A similar critical theory 
for this TQCP consisting of a single species of gapless 
Fermion can be calculated in _D = 2 in an analogous 
manner. In Z? = 1 this is the same as the critical the- 
ory of the Ising model in a transverse field, which in one 
dimension can be mapped into Eq. (fTO|) by a non-local 
Jordan- Wigner transformation [1]. 

Note that the nodal quasiparticles constitute a quan- 
tum critical phase which is essentially a non-interacting 
gas of two species of chiral Majorana fermions defined 
by 7|(fc) = 7i(— fc),i = 1,2. The gaussian critical point 
implied in Eq. (fTO]) is stable against interactions because 
all four-fcrmion interaction terms can be shown to be ir- 
relevant [1]. Furthermore, the effective action in Eq. (fTO|) 
involves only one species of regular fermion that be- 
comes gapless at the critical point. This is a key dif- 
ference between the Dirac spectrum found here and the 
analogous Dirac spectrum of the nodal quasiparticles in, 
say, dn^i_yi superconductors or HgTe quantum wells (27| 
where there are two species of gapless fermions corre- 
sponding to the spin degeneracy. Thus, our system avoids 
the fermion doubling theorem consequently giving rise 
to Majorana fermions and topological superconductivity 
whereas these other systems do not. 

In the vicinity of the TQCP, with the effective critical 
theory in Eq. (fTO|). the dissipative susceptibility xi^yi) 
for small a; and 8 (w, 5 much smaller than the gap at 
k = kp) takes a universal scaling form. This can be 
obtained by defining the rescaled variables, k = fc/T, 
= Ek/T, w = w/T and (5 5/T. The energy in 

Eq. [2] then takes the form E — \/ v'^lz^ + 5"^, while the 
matrix element fc|m, fc') p is invariant under the re- 
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FIG. 4. Dissipative susceptibility X2{T,u,5) shows scaling 
in the vicinity of the topological critical point 5 = F — Fc = 0. 
Appropriately scaled susceptibility data for ui/T = 0.5, cal- 
culated using Eq. ^ for different values of uj and T, coincide 
and collapse on the scaling function fEg. llip in the region near 
the critical point (5 ~ 0). The slight asymmetry of the data 
about 5 = is due to the asymmetry of the zero-temperature 
gap E(i{V) about F = Fc (5 = 0). 



scaling. The dissipative susceptibility in the vicinity of 
the TQCP takes the scaling form, 

X2^T^f{oj/T,6/T) (11) 

where the scaling function / is given by, 

m.n— ± * 



fi^,y)= E / ■ 



Snxv^k' 



-\{n,k\m,~kyCk + k')\'^ 



tanh I 



tanh 



E.., 



(12) 



Here E^j^ = ± 



fc^ -|- y2 and k' in the integrand is im- 



plicitly given by the equation E^ 

The existence of such a scaling function suggests that 
for a; < T {uj/T 0), X2/T^ = f{uj/T,6/T) depends 
only on d/T. In this limit, the scaling function in Eg. [T2l 
becomes f{0,y) ~ (l/7r)y^e~'', which becomes apprecia- 
ble only when the thermal energy scale ksT ~ S. This 
justifies our identification of /cbT ~ Eq curves (near the 
TQCP, Eo = r - Tc = 6, see Eq. O with the appro- 
priate crossover curves in the finite temperature phase 
diagram in Fig. [T]d. For to ^ T, the argument w/T 
in the scaling function / approaches 00 and must drop 
out. Thus, in this limit X^/T"^ is a function only of 
{ijj/T){5/T)-^ = uj/S. This implies that near T = 0, the 
dissipative susceptibility becomes appreciable only when 
Lo approaches 5 (and above). This is shown in Fig. [S] 
In general, the validity of such a scaling function rep- 
resentation becomes clear from Fig. IH where X2/T'^ as 
calculated from Eq. IH] is plotted for fixed w /T and com- 
pared with the scaling function Eq. [T^l The collapse of 
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the data for different uj and T (but with the same fixed 
ratio w/T) on the scaUng function (dashed curve) plotted 
as a function of 5/T is a clear and definitive experimen- 
tal signature of the topological quantum critical point at 

r = re. 

Finite-T crossover in supercurrent response: So 

far we have concentrated only on the dissipative part of 
the Josephson response of the nanowire bridge between 
the two superconducting leads in Fig. [51 The dissipative 
part allows us to access the nodal quasiparticles across 
the QC regime. The low-T finite-frequency response in- 
dicates an underlying gap collapse separating a fully- 
gapped s-wave superconductor at low T from another 
fully-gapped s-wave superconductor at high F that has 
identical broken symmetries (pair potential A remains 
the same). Although such a gap collapse at isolated 
points in the momentum space gives indications of an un- 
derlying TQCP, the evidence nonetheless is still circum- 
stantial. In particular, no information about the specific 
topological nature of the critical point and the high-F 
TS state can be derived from the behavior of the dissi- 
pative response across the QC regime. We now consider 
an indicator of the topological character of the underly- 
ing zero-temperature critical point from the behavior of 
the corresponding supercurrent response across the QC 
regime at finite temperatures. 

The experimental set up in Fig. [2] used to measure 
the finite-frequency dissipative response can be used to 
measure the supercurrent response as well. When the 
voltage V is time-dependent V{t) — y(cj)e"^*, the super- 
current response is given by the quantity Xi('^)- This is 
expected to behave as Lj/io throughout the phase di- 
agram {Lj is the effective inductance associated with 
the Josephson junction), and therefore cannot distin- 
guish between the NTS and TS states. Let us there- 
fore consider the limit uj <C V{uj) (where the linear re- 
sponse function Xi(<^) longer determines the super- 
current response) and in particular the case when V{t) 
is time-independent. Even if V{t) is time-independent 
the Josephson phase $ still linearly depends on time, 
$(t) = Vt. In this case, the Josephson current I{t) 
through the wire bridge connecting the superconducting 
leads in Fig. [2] should oscillate sinusoidally with a fre- 
quency determined by the applied DC voltage. Below we 
refer to this frequency as the AC Josephson frequency 
and derive its crossover behavior across the QC regime 
at low and finite T . At T = 0, such an experiment has 
been proposed U, 12, [il, [1^ to uniquely identify the 
topological character of the nanowire TS state. We show 
that the AC Josephson frequency also shows a crossover 
across the fc^T ^ Eq crossover line separating the QC 
regime and the TS state in Fig. [TJ At this crossover a 
peak at a fractional frequency in AC Josephson effect 
becomes dominant and the conventional Josephson fre- 
quency, characterizing the NTS state at small F, only 
makes a sub-dominant contribution. Tracking the AC 



Josephson response across the QC regime on the same 
set up as that for the dissipative response can uniquely 
identify the underlying TQCP and the subsequent TS 
state at large F. 

The current I(t) across a nanowire junction between 
two superconducting leads with a voltage difference V{t) 
is given by the relation I{t) = ^ , where Etot 
is the total energy in the system. The energy of a 
Josephson junction can be decomposed into two parts, 
Etot — E 



Eqp, where Ejunc is the energy stored 



in the localized states around the junction while Eqp is 
the energy dissipated into the quasiparticles that prop- 
agate away from the junction. The quasiparticle con- 
tribution to the current at low-frequency is dissipative 
and can be described by the dissipative response func- 
tion given earlier. The energy stored in the junction, 
Ejunc, is a function of both the energy of the Andreev 
bound state eABs{^), which depends on the phase dif- 
ference <&, and the occupation number n — 0,1 of the 
Andreev bound state. Specifically, shifting the occu- 
pation number n{t) of the junction switches the sign 
of Ejunc (because of particle-hole symmetry) so that 



Ejunc{'^,n) 



-l)"eABs(^')- The usual Josephson 



supercurrent carried by the Andreev state in the ground 
state (i.e. n = 0) is dissipationless and is given by 
Isc{^) = "^^^fp^i while the true current through the 
junction J($ , n) = 



a-Ej,.,ic(^.") 



of (-1)" [2c 
therefore 



has an additional factor 
For a fixed DC voltage, $(t) = Vt and 



/(t) =/($(*), n(i))-(-l)»(*)j5c(^t) 



(13) 



In the NTS state of Eq. ([T]) the numerical results for eABS 
can be fit by. 



eABs(*) = EoVTTDcc^, 



(14) 



while in the TS state it crosses zero energy and can be 
approximated by 



eABs($) 



r- * 

EqW D cos — , 



(15) 



where L* < 1 is the effective interface transparency [30|. 
While these spectra (together with the particle-hole sym- 
metric partners) look similar, particularly in the regime 
-D — > 1, they are fundamentally different in terms of 
fermion parity. In particular, for a fixed fermion num- 
ber n = in the NTS state the energy of the junction 
Ejunci^jn = 0) remains negative for all values of the 
phase while in the topological state the state with 
a fixed fermion parity n crosses zero energy. The AC 
Josephson current in the presence of a DC voltage quali- 
tatively distinguishes between the NTS and TS states. In 
the NTS state we get a ground-state (i.e. n = 0) current. 



I{n{t) = 0,i) = A 



D sin Vt 
^/l + Dcos Vt' 



(16) 
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FIG. 5. Spectral decomposition of current {|/(cj)p) as a func- 
tion of rescaled frequency uj/\S\ where 5 = F — Fc at several 
values of 5 characterized by 5/T. The applied DC voltage is 
taken to be V = 5/2 so that the AC Josephson current in the 
non-topological phase (i.e. 5 < Q) shows a peak a.t ui = 5 /2 
characteristic of the conventional AC Josephson effect. At 
and above 5/T ~ 1 (i.e. the crossover curve fesT ~ Eo for 
(5 > 0) a peak appears at uj = 5/ A which signals the ap- 
pearance of the fractional Josephson effect and an underlying 
T = TQCP at 5 = 0. 



which has harmonics only at multiples oiV. In contrast 
the current in the topological phase is given by the rela- 
tion 



I{n{t) = 0, t) = A\/i? sin ^ 



(17) 



and has a frequency oiV/2. 

The above picture for the supercurrent in terms of An- 
dreev bound states is only valid for biases that are smaller 
than the bulk gap of the system i.e. V < Eq. For higher 
biases, V ^ Eq, the harmonically varying superconduct- 
ing pairing term Ae*^* can excite a quasiparticle out of 
the Andreev bound state into the quasiparticle gap mak- 
ing the Andreev bound state ill-defined. Therefore the 
definition of the Josephson effect is valid for voltages sig- 
nificantly smaller than the gap. 

For voltages that are smaller than the gap (i.e. V < 
Eq), the Josephson current has a sign that is determined 
by the occupation of the Andreev bound state, namely 



(18) 



where n{t) is the occupation of the Andreev bound state. 
Here we assume n{t) to change instantaneously between 
the values and 1 and vice versa. Fermion parity conser- 
vation requires that such a change be accompanied by the 
emission of a quasiparticle which costs energy e^p > Eq. 
Thus the transition from the state n = to n = 1 requires 
an energy of Cqp + eABsi^) while the reverse transition 
requires an energy Cqp — CABsi^)- zero temperature 
such excitations are forbidden and the Josephson effect 
operates without fluctuations of n{t) as discussed before. 
At flnite T, the energy required for the above transitions 



is provided by thermal fluctuations either in the form of 
phonons or quasiparticles. In the following we refer to 
such excitations as phonons even though, the results will 
apply to more general excitations. For simplicity, we as- 
sume that a phonon above the required energy threshold 
will flip n(t) whenever it is energetically allowed. For 
temperatures T <^ Eq, the density of phonons with en- 
ergy e incident at the junction is given by the Bose- 
Einstein distribution nBE(f) = e'/'^-i • Assuming that 
the phonons move with the sound velocity, Vk, the rate 
of transitions at the junction for energy bigger than e is 
given by 



Pphonon 



de 



,V kie')Die') 



Te- 



rn 



for Eq > T. Here we have used that ^(e) = C-^T^ 

and Vk = ^^^iE^- expected the rate has the dimen- 
sions of energy. The corresponding flip rates can then be 
written as 

P{n{t) =0^1)= re-(-^°+^-^«^(*(*)»/'^ (20) 
P{n{t) = 1^0)= re-(-^°-^^«^(*(*)»/^. (21) 

In the low temperature limit {T < Eq) one can apply a 
voltage V such that Eo{1-VD) >V > Te'-^"/^. In this 
case the Josephson oscillation frequency is much bigger 
than the flip rate for n{t). For runs of I{t) where n{t) 
changes only after many periods of oscillation of Isc (t) , 
the Fourier transform should show a pronounced peak 
a.t uj — V for the regular Josephson effect (on the NTS 
side) and uj = V/2 for the fractional Josephson effect 
(on the TS side). In the temperature range T > Eq one 
is restricted to low voltages, V < Te~-^°^'^, since V is 
bounded by Eq . In this case the dynamics of the quasi- 
particle state n{t) is rapid compared to the phase and one 
can assume that n(t) is in local equilibrium, so that the 
resulting current I{t) will be a 27r-periodic function of the 
phase ^ = Vt. To show the crossover in the frequency 
dependence of the Josephson current across T ^ Eq, we 
calculate /(w) as the Fourier transform of I{t) in Eq. (14) 
and calculate the average (|/(w)p) where the average is 
taken with respect to random realizations of n(t) accord- 
ing to Eqs. (16,17). In Fig. [5] we show the evolution of 
the frequency dependence of the AC Josephson response 
at flxed low T for various values of S/T (a horizontal cut 
across the QC regime in Fig. [T|d). The flnite-T Joseph- 
son response shows a pronounced crossover across the 
ksT ^ Eq curve separating the QC regime from the TS 
state at large F. At this crossover, frequency peaks in- 
dicating an underlying zero-T Josephson period-doubling 
transition starts making a dominant contribution and the 
usual Josephson frequency becomes sub-dominant. Such 
a crossover in the superfluid response can serve as an un- 
ambiguous marker of an underlying zero-T TQCP in the 
semiconductor heterostructure. 
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Summary and Conclusion: TQCPs separate two 
macroscopic ground states which have the same sym- 
metries and hence cannot be distinguished by an order 
parameter. Consequently, the topological and the non- 
topological states on the two sides of a TQCP cannot in 
general be distinguished by any bulk measurement, e.g., 
no thermodynamic quantity diverges at the TQCP. To 
solve this problem, we propose the use of specific aspects 
of the topological phase transition itself as an identifier. 
We consider the TQCP in a spin-orbit coupled semicon- 
ductor {e.g., InAs) thin film or nanowire on which s-wave 
superconductivity is proximity induced. The TQCP in 
this case is tuned by an external Zeeman splitting F, and 
for large F the ground state of the system is a topo- 
logical superconductor. We ask if straightforward bulk 
transport measurements can help identify the emergence 
of the TS state with increasing values of F. We show 
that this is indeed possible and establish that the finitc- 
T AC Josephson impedance along with the supercurrent 
response of the semiconductor can access the topological 
quantum critical regime which precedes the TS state in 
the finite-T phase diagram (Fig. [T|). 

We do this by first identifying the entire QC regime 
with a bulk s-wave superconductor coexisting with gap- 
less nodal fermions at fc = 0. Since both the non- 
topological superconducting state at low F and the TS 
state at high F are fully gapped, both these states are de- 
void of quasiparticles at sufficiently low T,uj <^ Eq, where 
Eq{T) is the Zeeman-tunable single-particle energy gap. 
The real part of the inverse AC Josephson impedance 
(X2)j which gives the dissipative response of a nanowire 
contacted by two superconducting leads (Fig. [2]), is there- 
fore negligible away from the TQCP for cj,T less than 
Eq in both the NTS and TS states. As the TQCP is ap- 
proached with increasing F, the Josephson impedance for 
a given frequency picks up as the frequency becomes com- 
parable with the decreasing values of the energy gap. In 
the low temperature QC regime in the finite-T phase dia- 
gram, the impedance follows a scaling function involving 
T,uj, and (5 = F — Fc, which can be experimentally veri- 
fied to help reveal the underlying quantum critical point. 
When the Josephson impedance for a given frequency de- 
creases again with increasing values of F past Fc, it indi- 
cates the emergence of the TS state on the high-Zeeman- 
field side of the critical point. The finite-T superfluid 
response in the set up of Fig. [2] also shows a pronounced 
crossover at the ksT = Eq curve on the large-F side of 
the phase-diagram (Fig.[TlD). At this crossover, the usual 
Josephson current frequency in response to a DC voltage 
becomes sub-dominant and pronounced peaks at a frac- 
tional frequency dominates the spectral decomposition of 
the Josephson current (Fig. [5]). Such a finite-T crossover 
in the superfluid response can be observed in the same 
set up as in Fig. [2] and can serve as an unambiguous 
indicator of the underlying T — TQCP in the semicon- 
ductor heterostructure. Our work demonstrates that a 



bulk experimental characterization of TQCP may indeed 
be possible although it may require a careful analysis of 
specific properties of the topological phase transition as 
we have carried out here for the TS state in semiconduc- 
tor heterostructure systems with the real significance of 
our work lying in the fact that no TQCP has yet been 
clearly identified experimentally in any system. 
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